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Abstract. This is a rather personal review of the problem of self-avoiding walks and polygons. After defining 
the problem, and outlining what is known rigorously and what is merely conjectured, I highlight the major 
outstanding problems. I then give several applications in which I have been involved. These include a study of 
surface adsorption of polymers, counting possible paths in a telecommunication network, hitting probabilities of 
SAWs in a rectangle, and the modelling of biological experiments on polymers. I hope to show that SAWs are not 
only of intrinsic mathematical interest, but also have many interesting and useful applications. 



1. Introduction 



The problem of self-avoiding walks is one of deceptive simplicity of definition, hiding malevolent difficulty of 
solution. The problem was introduced by two theoretical chemists, Orr [35] and Flory |TS] , as a model of a 
polymer in dilute solution. It soon became an interesting combinatorial model to mathematicians, and a canonical 
model of phase transitions, of interest to mathematical physicists. It is also a simple model of a non-Markovian 
process. Attempts to count the number of SAW have led to the development of new algorithms, with widespread 
applicability, while many more applications were discovered. These include application to the design of telephone 
networks, the folding and knotting of biological molecules, and a variety of chemical phenomena. Attempts at a 
solution have driven several mathematical advances, including developments in stochastic differential equations 
and probability theory. 

Nearly 70 years after the model was proposed, we have a huge amount of numerical information, a substantial 
amount of exact information - that is to say, results that are universally believed, but remain unproved - and 
a very small body of rigorous results. In contrast, some other canonical models of phase transitions, such as 
the Ising model, the Potts model and percolation have either been solved (the Ising model) or much has been 
rigorously proved. In this short article I will outline the development of the subject, give some applications, and 
show that we appear to be on the verge of some major breakthroughs, which will result in proofs of much of the 
exact, but unproved, information that currently exists. Unfortunately all the exact and conjectural information 
we have applies only to the model on a two-dimensional lattice. In the case of three dimensions, we only have 
numerical results. Except where otherwise stated, this article will discuss the two-dimensional situation. 

Finally, a veritable treasure trove of rigorous results would be unlocked if we could prove that, in the large size 
limit, more precisely the scaling limit, that two-dimensional random SAWs are described by one of the SLE 
processes (Schramm-Loewner Evolution), which, in the past 20 years, have been proved to describe the limit of 
several discrete models in combinatorics and statistical physics. Indeed, Lawler, Schramm and Werner proved 
that, if the scaling limit of SAWs exists and is conformally invariant, then this limit has to be SLE 8 / 3 . This 
has been checked via simulations by Kennedy |25j . This would in particular imply the conjectured values of the 
exponents g and v, which are defined below, in the two-dimensional case. We explain this in more detail in the 
conclusion. 
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2. What is known and what isn't. 

2.1. Self-avoiding walks 

A self-avoiding walk (SAW) of length n on a periodic graph or lattice £ is a sequence of distinct vertices 
Wo, iui, . . . , w n in £ such that each vertex is a nearest neighbour of its predecessor. In Figure [T] a short SAW on 
the square lattice is shown, while in Figure [2] a rather long walk of 2 25 steps is shown (generated by a Monte Carlo 
algorithm (10115]). 



FIGURE 1. A self-avoiding walks on the square lattice. 




Figure 2. A typical two-dimensional SAW of 2 25 steps on the square lattice - courtesy of Nathan Clisby. 



2.1.1. How many self-avoiding walks are there? Two obvious questions one might ask are (i) how many 
SAWs are there of length n, (typically defined up to translations) denoted c n , and (ii) how big is a typical n-step 
SAW? Indeed, how might we measure size? A third important, but less obvious question, asks "what is the scaling 
limit of SAWs?" 

Frequently one rather considers the associated generating function 

c(x) = c ^ n - 

n>Q 

To see the difficulty of this problem, the reader is invited to try and calculate the first few terms c„ on Z 2 . We 
take Co to be 1, then c\ = 4 as a one step walk can be in any of 4 directions. Then C2 = 12, C3 = 36 and C4 = 100. 
It is at the stage of 4-step SAWs that the self-avoiding constraint first manifests itself, and the problem becomes 
increasingly difficult thereafter. 

As we prove below, c n grows exponentially. Accordingly, an enormous amount of effort has been expended over 
the last 50 years in developing efficient methods for counting SAW. For the square lattice, Jensen [?] has extended 
the known series to 79 step walks, for which he finds c 79 = 10194710293557466193787900071923676. Methods for 
calculating these astonishing numbers are quite complicated (see |16| , Chapter 7) , but the best current algorithm 
still involves a counting problem of exponential complexity, of about 1.3™ (while a direct counting algorithm 
would have complexity 2.64"). 



SELF-AVOIDING WALKS 



3 



One of the few properties one can readily prove, by virtue of the obvious sub-multiplicative inequality c n+m < 
c„c m , is that the number c n grows exponentially. From this inequality it follows that 

/i:= lim cl/ n = inf cl/ n 

n— foo n 

exists [33 , and further that c n > /i™. 



However even the value of this "growth constant" [i is difficult to calculate exactly. Only in 2010 was \x for 
one two-dimensional lattice, the honeycomb lattice, actually proved by Duminil-Copin and Smirnov [12 to be 

see Sectional. For other lattices in two dimensions, and all lattices in higher dimensions, we only have 
numerical estimates. For example, for the square lattice the best current estimate is [i = 2.6381585303±2 x 10 -11 , 
given by Clisby and Jensen [23 . 



In fact it is believed that, for dimensionality d > 1 and d 7^ 4, 

c n ~ const. x/x n n 9 . 

The critical exponent g is believed to depend on the dimension, but not on the details of the lattice. In particular, 
it is predicted to be a rational number, namely 11/32, in two dimensions. In three dimensions, the best estimate 
we have is g = 0.156957±0. 000009 given by Clisby |11J . There is no reason to believe that this number is rational. 



Despite these accurate estimates, we still cannot even prove the existence of this exponent for d < 5, let alone 
establish its value rigorously. For d > 4 the higher dimensionality means that the self-avoiding restriction is 
less confining than in lower dimensions, and indeed has no effect on the dominant asymptotic behaviour, with 
the result that the SAW behaves as a random walk. More precisely, Hara and Slade [211 |2"U] have proved that 
g = in this case, and that the scaling limit is Brownian motion. In four dimensions the above expression for c n 
must be modified by an additional multiplicative factor (logn) 1//4 , with g = 0. The appropriately rescaled walk 
is also expected to have Brownian motion as its scaling limit. These assertions for the four-dimensional case are 
believed to be true, but no proof exists. See [7] for a discussion of this case. Bounds established 50 years ago 
by Hammcrslcy and Welsh [19 have hardly been improved upon. They proved that, for SAW in dimensionality 
d > 2, 

The lower bound follows immediately from sub-additivity, while the upper bound depends on an unfolding of the 
walk. The number of possible unfoldings can be bounded by the number of partitions of the integer n, which 
has the exponential behaviour given above. Note that the existence of a critical exponent would imply behaviour 
/i n e Klog ", which is rather far from the upper bound. A year later, Kesten (26] slightly improved the upper bound 
for d > 2 to 

< „n Kr> 2 /< d + 2 >logn 



2.1.2. How large is a typical self-avoiding walk? Another important measure of SAW is the average size 
of a SAW of length n, taken uniformly at random. The most common measure is the mean-square end-to-end 
distance, which is believed to behave as 

E„(|w„| 2 ) ~ const, n 2u , 

(for lattices in dimensions other than 4), where v is another critical exponent. Again, its existence hasn't been 
proved for d < 5, but it is accepted that for two-dimensional lattices v = 3/4. In three dimensions the best 
numerical estimate is v — 0.587597 ± 0.000007 [9 . In four dimensions it is believed that v = 1/2, and again one 
expects a multiplicative factor (logn) 1 / 4 . Finally, for d > 4 it has been proved [2T] that v = 1/2. Rigorous results 
about E n (|«; n | 2 ) are almost non-existent. It would seem intuitively obvious that 

cn < E n (\w n \ 2 ) < Cn 2 -\ 

but only this year, in an important calculation, has substantial progress in establishing the upper bound been 
made by Duminil-Copin and Hammond [13J for SAW in dimension two or greater. They prove that walks are 
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sub-ballistic, in the sense that there is an exponentially small probability of them having any given positive speed. 
But technically, the upper bound still remains an open probleirj^] 



2.2. Self-avoiding polygons 



If the end-point of a SAW is adjacent to the origin, an additional step joining the end-point to the origin will 
produce a self-avoiding circuit. If we ignore knowledge of the origin, and distinguish circuits only by their shape, 
we refer to self-avoiding polygons (SAP). On the square lattice, the first non-zero embedding of a SAP is the unit 
square, of perimeter 4 and area 1. There are two 6-sided polygons of area 2, and seven 8-sided polygons, shown 
in Figure [3j one of which has area 4, and six of which have area 3. 



FIGURE 3. All seven 8-sided polygons on the square lattice. 

Clearly, SAPs are a subset of SAWs. They are a particularly interesting subset for at least two reasons. Because 
the conjectured exponents for SAPs (discussed below) are integers or half-integers (which is not the case for 
SAWs), it is hoped that this means the underlying solution for the SAP case is simpler. Secondly, by including a 
second parameter, that of area, SAPs can be used to model a range of biological phenomena, such as cell inflation 
and collapse |14| . 

Denote by p m the number of SAPs of perimeter m, by a n the number of SAPs of area n, and by p m . n the number 
of SAPs of perimeter m and area n. We can define two single variable generating functions, for perimeter and 
are80 respectively, and a two- variable generating function, as follows: 



P(x) = 5> m x r ' 

m 

A(q) = ^2 a n q n 



V{x,q) = Y J 



n r m a n 



Hammersley [TB] proved that the number of SAPs, like SAWs, grows exponentially; more precisely 

n . l/2m 

fx = km p 2 ' m . 

m— foo 

While it is far from obvious, Hammersley also proved that the growth constants [i that arise in the polygon case 
and the walk case are identical. While unproved, a much stronger result is widely believed to hold, namely that 

p m ~ const. X(i m m a ~ 3 (1) 

where a is a critical exponent ^ The exponent a is related to the exponent v defined above through the hyper- 
scaling relation dv = 2 — a. This equation has not been proved, but follows from physical arguments, and of 

4 am grateful to Alan Hammond for clarification as to precisely what has been proved. 
2 Clearly the area of a polygon is a concept peculiar to the two-dimensional case. 

^Note that p2m+i = for SAP on Z d , as only polygons with even perimeter can exist on those lattices. For such lattices the above 
asymptotic form is of course only expected to hold for even values of m. For so called close-packed lattices, such as the triangular or 
face-centred cubic lattices, polygons of all perimeters greater than two are embeddable, so eqn. |TJ stands as stated. 
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course the assumption that the exponents exist. It therefore follows from the result for v quoted above that in 
three dimensions a = 0.237209 ± 0.000021. 

For polygons there is a second growth constant, and exponent, associated with the area generating function. By 
concatenation arguments it can be readily proved that 

A = lim a l J n 

n— voo 

exists. It is also generally accepted, but not proved, that 

a n ~ const x X n n T . 

Unfortunately we only have numerical estimates of A and r [16] . However for two-dimensional lattices r is believed 
to be —1, corresponding to a logarithmic singularity of the generating function. That is to say, 

A(q) ~ const x log(l — Xq), 

so that a n ~ const x X n /n. 

Of great interest is the two- variable generating function V(x,y). From this, we can define the free energy 



n(q) = lim — log 



It has been proved |14| that the free energy exists, is finite, log-convex and continuous for < q < 1. For 
q > 1 it is infinite. The radius of convergence of V(x, q), which we denote x c (q), is related to the free energy by 
x c{q) = e~ K ( q ' . This is zero for fixed q > 1. A plot of x c (q) in the x — q plane is shown (qualitatively) below. For 
< q < 1, the line x c (q) is believed to be a line of logarithmic singularities of the generating function V(x,q). 
The line q — 1, for < x < x c (l) is believed to be a line of finite essential singularities [14]. At the point (x c , 1) 
we have more complicated behaviour, and this point is called a tricritical point. 




Around the point (x c , 1) we expect tricritical behaviour, so that 

V {sma \x,q)~ (l~q) e F{{x c -x)(l- q r*) as (x,q)^(x c ,l-). (2) 

Here the superscript (sing) means the singular part. There is an additional, additive part that is regular in the 
neighbourhood of (x c , 1). 



For self-avoiding polygons, in a series of papers, Richard and co-authors |37|, [38j [36] have provided abundant 
evidence (but no proof) for the surprisingly strong conjecture that 

1 



F(s) 



2?r logAil --(4Ao)s S ) +C(q). 
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Here Ai(x) is the Airy function, and C(q) is a function, independent of x, that arises as a constant of integration. 
Then equation ^ implies that the exponents introduced in equation Q are = 2/3 and 9 = 1. Here Aq is a 
constant known only numerically. 

In the next section we give the proof due to Duminil-Copin and Smirnov of the exact growth constant for the 
honeycomb lattice. In the following section we give three examples of applications of SAWs to other areas of 
science, and in the conclusion we give more detail of recent developments that we hope point the way to future 
breakthroughs. 



3. The honeycomb lattice 

As mentioned above, a breakthrough was achieved in 2010, when Duminil-Copin and Smirnov [T^j proved that 

the growth constant on the honeycomb lattice is fi = v2~+V2, as predicted by Nienhuis [53], using compelling 
physical arguments from conformal field theory, 30 years previously. The argument is, in hindsight, so simple, 
and the result so important, that we sketch it here. 




Rh,l 



I 



Figure 4. A trapezoid T on the honeycomb lattice. 



We consider SAWs that start from a point a located on the left side of a trapezoid T of width I and height h, as 
shown in Figure [4j For p a mid-edge of T, let F(p) be the generating function of SAW w that end at p, weighted 
by the number of vertices v(w) and the number of turns T(w) (a left turn counts +1, a right turn —1): 

F(p) = F(p;x,a):= J] i vM e« aT ('">. 

For instance, the walk of Figure [4] visits 17 vertices, makes 10 left turns and 7 right turns, so that its contribution 
to F(p) is x 17 e 3la . Then, if v is any vertex of T and pi,P2,P3 are the three mid-edges adjacent to it, the following 
local identity holds: 

(pi - v)F{ Pl ) + (pa - v)F(p 2 ) + (pa - v)F(p 3 ) = 0, (4) 

provided x = x c := l/y2~+ \/2, which is the reciprocal of the conjectured growth constant, and a = — 57r/24. 
(We consider that the honeycomb lattice is embedded in the complex plane C, so that pi — v is a complex number). 
This identity is easily proved by grouping as pairs or triplets the SAWs that contribute to its left-hand side, as 
depicted in Figure [5] One then checks that the contribution of each group is zero. 
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Figure 5. A very local proof of the local identity Q. 

If we now sum over all vertices v of T, then due to the terms (pi — v), all terms F(p) such that p is not a 
mid-edge of the border disappear. After a few more reductions based on symmetries, one is left with 

^cos L h: e(x c ) + M h j(x c ) + R h ,e{x c ) = 1, 

where Lh^(x) (resp. Rh^(x), Mh^{x)) is the generating function of SAWs w that end on the left side (resp. right 
side, top or bottom) of T, weighted by the number of vertices v(w). 

By letting h and then I tend to infinity, Duminil-Copin and Smirnov derived from this identity that the generating 
function of SAWs diverges at x c , but converges when x < x c . This means that its radius of convergence is x c , so 
that the growth constant is l/x c = y2 + y/2. 

Unfortunately these ideas do not generalise to SAW on the square or triangular lattices, for which we only have 
accurate numerical estimates for the growth constant fj,. 

4. Applications 

One reason that SAWs and SAPs are so extensively studied, apart from their intrinsic mathematical interest, is 
that they model many problems that arise in other fields. The first such example we will consider extends the 
proof given above to the situation where the SAW can interact with a surface. The second example considers 
SAWs crossing a square, with application to telecommunication networks, and the third example models some 
recent biological experiments where strands of DNA (a polymer) are pulled from a wall with optical tweezers. 

4.1. Walks attached to a surface 

The interaction of polymers with a surface is scientifically and industrially an important phenomenon. A common 
example is the adherence of paint to a surface, clearly an industrial process of considerable significance. To model 
such phenomena requires the inclusion of an interaction term between the polymer and the surface. To achieve 
this, we add a weight y to vertices in the surface, as shown in Figure [fj] below. In physics terms, y = e _e / fc£iT 
where e is the energy associated with a surface vertex, T is the absolute temperature and ks is Boltzmann's 
constant. It is known that the growth constant fi — l/x c for such walks is the same as for the bulk case. 

Let c+(i) be the number of half-plane walks of n-steps, with i monomers in the surface, and define the partition 
function (or generating function) as 

n 

C+(y) = J24(*W- 

If y is large, the polymer adsorbs onto the surface, while if y is small, the walk is repelled by the surface. 
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Proposition 1. For y > 0, 

l*(y) := lim C+(y) 1/n 

n—too 

exists and is finite. It is a log-convex, non- decreasing function of logy, and therefore continuous and almost 
everywhere differentiable. 

For < y < 1, 

M(y) = = 

Moreover, for any y > 0, 

H{y) > max(/J,,y/y). 

This behaviour implies the existence of a critical value y c , with 1 < y c < [i 2 , which delineates the transition from 
the desorbed phase to the adsorbed phase: 



= a 1 -ify<y c , 

>M if V> Ve- 



in 1995 Batchelor and Yung pQ extended Nienhuis's [33] work to the adsorption problem just described, and 
making similar assumptions to Nienhuis conjectured the value of the critical surface fugacity for the honeycomb 
lattice SAW model, to be y c = 1 + y/2. In 2012 this was proved by Beaton, Bousquet-Melou, de Gier, Duminil- 
Copin and Guttmann (3j, and here we will sketch their proof. 

Take the same trapezoid as above, now called «St;.l, and add weights to the vertices on the /3 boundary, as 
shown in bold in the figure below: Then we find the corresponding identity between generating functions, with 




FIGURE 6. Trapezoidal domain St,l with vertices on the right-hand side wall, shown in bold, 
carrying a weight y. 

y* = 1 + V2, to be 

1 = cos ( — J A T ,L(x c ,y) +cos (^-J E T , L {x c ,y) + _ ^ B T ,L(x c ,y) 

where y is conjugate to the number of visits to the /3 boundary. The generating functions At^l, Bt.l and Et,l 
are two- variable generalisations of those defined in the previous section. To prove the conjecture we need to show 
that y c = y*. 

It is safe to take L —¥ oo so that the trapezoid becomes a strip. The identity then becomes 



(t) At ^ Xc ' V) + C ° S (l) Et ^ Xc ' y) + yly* V \) Bt ^ y) - 
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It is then straightforward to prove that 

(i) E T {x Cl y) = 0ior0<y<y*, 

(ii) y c >y*, 

(hi) lim.r_j.oo At{x Ci d) = A(x c ,y) = A(x c ) is constant for < y < y* . 
If we now write 

cos(37r/8)A(a; c ,y) = 1-5 

then the above identity reduces to 

B(x c ,y) = lim B T (x c ,y) = — 

t->oo y y 

and in particular 

B{x c ,l) =6. 

Proposition 1. If S = then y c = y* . 

The proof uses a decomposition of A walks in a strip of width T into B walks in that same strip, and gives rise 
to an inequality. In particular, for y < y c — lim^-foo 2/T, 

n s , 1 y* -y 

< ax c 



B T (x c ,l)y(y*-lY 

If Bt(x c ,1) tends to 0, this forces y* > y c , otherwise the right-hand side would become arbitrarily large in 
modulus and negative as T — » oo for y* < y < y c . 

Together with y* < y c , this establishes y c = y* = 1 + y/2 and completes the proof of the proposition. 

The proof that S = is complicated, and unlike most other proofs we have given is almost totally probabilistic. 
It is unrealistic to give any details, but in essence one first uses renewal theory to show that <5 _1 is the expected 
height of an irreducible bridge, which is a SAW that crosses the strip from left to right, and cannot be expressed as 
the concatenation of two or more smaller such bridges. Next one shows that, for irreducible bridges, E[width] < oo 
implies that K[height] < oo. Finally one shows that the assumption that E[height] < oo leads to a contradiction, 
from which the desired result that <5 = readily follows. 



4.2. Walks crossing a square 



Some years ago I was asked by a telecommunications engineer to help him with the following problem: His 
company had a square grid of nodes, connected by wires, and phone-calls could be routed from the bottom 
left-hand corner to the top right-hand corner of the grid. He wished to know how many such routes there were, 
as this determined the carrying capacity of the network. 

After some discussion we agreed that this was simply the question how many distinct SAWs are there on a square 
grid of side-length L originating at (0,0) and ending at (L,i)? The problem as stated was first considered by 
Knuth in 1976, who gave a Monte Carlo estimate for the number of paths for L = 10, a result we now know 
exactly. The problem was generalised by Whittington and the author (33] to include a weight x associated with 
each step of the walk. This gives rise to a canonical model of a phase transition. For x < l//i the average length 
of a SAW grows as L, while for x > it grows as L 2 . Here fj, is the growth constant of unconstrained SAWs 
on the square lattice, defined above. For x — l//i numerical evidence, but no proof, was given that the average 
walk length grows as L 4 / 3 . Let c n (L) denote the number of walks of length n. Clearly c n (L) = for n < 2L. We 
denote the generating function by Cl(x) := ^ n c n (L)x n . The answer to the original question is ^ n c n (L). 
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Subsequently, Madras [32 proved a number of relevant results. In fact, most of Madras's results were proved 
for the more general d-dimensional hypercubic lattice, but here we will quote them in the more restricted two- 
dimensional setting. 

Theorem 2. The following limits, 

Hl(x) := lim C L (x) 1/L and fjt^x) := lim C L (x) llL \ 

are well-defined in K U {+oo}. 

More precisely, 

(i) fJ-i(x) is finite for < x < 1/fi, and is infinite for x > 1/fx. Moreover, < fJ-i(x) < 1 for < x < 
and Hx(l/fj,) = 1. 

(ii) (Jq{x) is finite for all x > 0. Moreover, fJ<2(x) = 1 for < x < and fJ,2(x) > 1 for x > 

In |43| the existence of the limit [iiix) was proved, and in addition upper and lower bounds on [ii{x) were 
established. 

The average length of a (weighted) walk is defined to be 

(n(x, L)) := J2 nc n {L)x n / £ c n (L)x n . (5) 

n n 

We define Q(x) as follows: Let a(x) and b{x) be two functions of some variable x. We write that a(x) = &(b(x)) 
as x — > Xq if there exist two positive constants Ki and k 2 such that, for x sufficiently close to Xq, 

k\ b{x) < a(x) < k-2 b(x). 

Theorem 3. For < x < 1/fJ., we have that (n(x,L)} = O(L) as L — > oo, while for x > 1/fi, we have 
(n(x,L))=e(Li). 

In |43| it was proved that (n(l,L)) = <d(L 2 ). The situation at x = is unknown. In we gave compelling 
numerical evidence that in fact (n(l//i, L)) = 6(i 1 / 1/ ) , where v = 3/4, in accordance with an intuitive suggestion 
of Madras in [32]. 

Theorem 4. For x > 0, define fi(x) = log/xi(a;) and f2(x) — \ogfJ,2(x). 

(i) The function f\ is a strictly increasing, negative-valued convex function of logx for < x < 
(ii) The function / 2 is a strictly increasing, convex function of logx for x > 1/fi, and satisfies < f2{x) < 
log (j, + logx. 

Some, but not all of the above results were previously proved in |43| . but these three theorems elegantly capture 
all that is rigorously known. In [B] an extensive numerical study was described, including exact enumerations 
up to squares of side 19. For the largest square there are exactly 1 523 344 971 704 879 993 080 742 810 319 
229 690 899 454 255 323 294 555 776 029 866 737 355 060 592 877 569 255 844 distinct paths! The number of 
such paths, as we have seen, grows as \ L ~ . In |6] it was also proved that 1.628 < A < 1.782 and estimated that 
A = 1.744550 ± 0.000005. 

4.3. Walks in a rectangle 

In 2002 L N Trefethen [IT] presented ten problems used in teaching numerical analysis at Oxford University. 
Problem number 10 asks: A particle at the centre of a 10 X 1 rectangle undergoes Brownian motion (i.e., 2-D 
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random walk with infinitesimal step lengths) till it hits the boundary. What is the probability that it hits at one 
of the ends rather than at one of the sides? 

Remarkably, it turns out [5] that a closed form solution to this question can be found, via the theory of elliptic 
functions and invoking some results of Ramanujan on singular moduli. The result is 

p e = ^ arcsin ((3 - 2V2) 2 {2 + Vb f (VlO - 3) 2 (5 1/4 - v^) 4 ) = 0.000000383758797925 .... 

In |17| . the question as to the corresponding problem where Brownian motion is replaced by the scaling limit of 
a self-avoiding walk (SAW) is treated. In fact, they addressed the equivalent question, which is the ratio of the 
probability that the walk hits the end (j> e ) to the probability that it hits a side. This is just 

Pe 
1 - Pe 

which has the numerical value 0.00000038375894519599411176841999126970034234598936 . . . for the random walk 
case. 

By using conformal mappings, it was shown in [T7| how to solve this problem for both Brownian motion and for 
the scaling limit of SAW, assuming that this is describable by SLE s / 3 , as discussed in Sections 5.1 and 5.2 below. 

Let D be a bounded, simply connected domain in the complex plane containing 0. We are interested in paths 
in D starting at and ending on the boundary of the domain. For Brownian motion, the distribution of the 
end-point is harmonic measure. Let ho{z) denote the density with respect to arc length (often called the Poisson 
kernel). If / is a conformal map on D that fixes the origin, and such that the boundary of f(D) is also piecewise 
smooth, then the conformal invariance of Brownian motion implies that the density for harmonic measure on the 
boundary of f(D) is related to the boundary of D by 

h D {z) = \f{z)\h f{D) {f{z)). (6) 

For SAW, Lawler, Schramm and Werner [31 predicted that the corresponding density of the probability measure 
p(z) transforms under conformal maps as 

p D (z) = c\f'(z)\ b p nD) (f(z)), (7) 

where 6 = 5/8 and the constant c is required to ensure that pd(z) is a probability density. Note that b is related 
to k by 6 = 3/k — 1/2, so 6 = 5/8 corresponds to k — 8/3 . If one starts the random walk or the SAW at the center 
of a disc, then the hitting density on the circle will be uniform. So the above equations determine the hitting 
density for any simply connected domain. 

The map / between the upper half-plane and a rectangle (see figure below), is given by a Schwarz-Christoffel 
transformation. For a > 1, let 

t(z)= r d( 

Jo v^eV^e' 

The rectangle has one edge along the real axis and is a midpoint of this side. So the corners can be written as 
±a/2 and ic ± a/2 where a, c > are the length of the horizontal and vertical edges, respectively. We have 



/(l) = a/2, /(-l) = -a/2, f(a) = a/2 + ic, /(-a) = -a/2 + ic, /(0) = 0. 
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We note that 



dx 



dx 



l \/ x 2 — ly/a 2 — x 2 



a = -K — , c = -K , a > 1. 

a \ a J a \ a ] 



Here K(x) is the complete elliptic integral of the first kind. By dilation invariance we only need concern ourselves 
with the aspect ratio a/c. So given an aspect ratio r, we have to find a such that 



J^{1 - x 2 )- 1 / 2 ^ 2 - x 2 y l l 2 dx 2K(£ 



J"(x 2 -l)- 1 /2( a 2 - X 2 )-V 2 dx k ( VaEE 



Setting k = 1/a, then c can be written K(fc') = K'(fe). Therefore 

c _ 1 _ K(l/q') _ K'(l/a) 
a~r~ 2K(l/a) ~ 2K(l/a)" 

Now the nome g is definecQby g = exp(— 7rK'/K), so 

2K(l/o) 



K'(l/a) 



and 



Then we have 



- g. 



J (9) 



so 

V " d 2 (e-a*/r)' 

where #j(g) = 9j(Q,q) is the Jacobi theta function. Evaluating this with any computer algebraic package gives 
the required value of a for any r > 1 instantly. 



Amramowitz and Stegun, p591 formula 17.3.17 
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Alternatively, we can achieve good accuracy by expanding the ratio of the above integrals around a = 1. This 
gives 

r = - (a\o S (2V2) - 21og(o - 1) + (a - 1) - |(q - l) 2 + A(a - 1)3 + (a - l) 4 

TT \ 8 24 

Solving this numerically for r = 10 gives a — 1.00000120561454706472212 .... Reverting this equation, retaining 
the first two terms gives 

a = 1 + %e^ rl2 + 32e- 7rr + 0( e - 37rr/2 ), (8) 
which, for r = 10 gives 19 significant digits. 

To proceed further, we need the preimage of the center of the rectangle, which by symmetry is on the imaginary 
axis, so write it as id. Either by symmetry arguments, or from the solution of the mapping equation, we readily 
find [17J that d 2 = a. 

For the random walk and the SAW in the half plane starting at id, it follows from ^ and (JgJ) that the (unnor- 
malized) hitting density along the real axis is (x 2 + d 2 )~ b = (x 2 + a)~ b . So letting pa denote the hitting density 
for a rectangle with the walk starting at the centre, we have by ^ and (Jo]) , 

{x 2 +a)- b ^\f{x)\ h p R {f{z)). (9) 

We require the ratio of the integral of p(z) along a vertical edge to the integral along a horizontal edge, 

Jo PR(a/2 + iy)dy 
Jl / a /2 Pnix)dx 

By a change of variable, setting u = f^ 1 (x) in the denominator and u = f~ 1 (a/2 + iy) in the numerator, this 
becomes 

j«p R {f{u))\f>{u)\du 

S\pR{f(u))\f>(u)\du 

Since f'(u) = (1 — u 2 )^ 1 / 2 (a 2 — u 2 ) -1 / 2 , the ratio of probabilities of a first hit on the vertical side to a first hit 
on the horizontal side, R(a, b) is 

C(u 2 + a)-Hu 2 - \)^/ 2 {a 2 - u 2 V b -V/ 2 du 

R(a, b) = ^ '— X V ; . 10 

y_ x {u 2 + a)- b (l - u^-Wiafl - u 2 )( b - 1 )/ 2 du 

Note that this result should hold for b 6 (1/4, 1], so applies more broadly than just to SAWs. 



The Brownian motion case corresponds to b = 1. The integrals in (10) greatly simplify, giving 



arctanfy^) — arctan(l/-v/a) 

RiaA) = ; — , ' . 

v ; 2arctan(l/Va) 

It is straightforward to calculate the asymptotic expansion of both the numerator and denominator, and hence 
their ratio. In this way we find 

R ^ i) = e -*r/2 + ^e-" - + O ( e- 3 ^/ 2 ) . 
For r — 10 this evaluates to 3.8375894519594. . . x 10~ 7 , which is correct to 13 significant digits. 

Unfortunately for 6^1 one cannot evaluate the integrals exactly, but one can evaluate the relevant integrals 
numerically. For SAWs b = 5/8 and one finds £(10,5/8) « 6.682989935 x 10~ 5 , some 200 times larger than the 
corresponding result for SAWs. 
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As for the random walk case, asymptotics also provides us with reasonable accuracy. Referring to equation (10 1, 
and noting that a is very close to 1 when the aspect ratio is 10, the denominator integral can be accurately 
approximated by 

i 21 U + 2) 



The numerator integral has a very narrow range of integration, and so can be approximated by 

2 Ji 

If we set u = 1 + t(a — 1), this becomes 

\{a - lfj\t{l - tr-^dt = 2 "V(f + r |f ) ( Q - 

Combining the above results, and using the result derived above, a — 1 « 8e~ 7rr / 2 , one finds, asymptotically 

2 2b T (1 , fA 2 

&) ~ V l2 M 2 L ' itr/2 . (11) 

r(i + |)r(|) 



For SAW equation (111 reduces to 

R(r,5/8) fa 1.2263431442e- 57rr/16 . 
For r = 10 this gives i?(10, 5/8) s» 6.6824528 x 10 -5 , which is accurate to 4 significant digits. 

A rather elaborate calculation using Mellin transforms gives the next term in the asymptotic expansion of the 
aspect ratio as 



R(r,b) 



for < 6 < 1, where A 



126+1 



A 



-birr/2 



K2 b+1 
bsm(^) 



-bur 12 



+ 4(6 - 1 + 2A)e-* r /' 2 + 0(e _6 * r ) 



r(l) 



Evaluating this for r = 10, b = |, we find 
E(10, 5/8) = 0.00006682989679 • • • 



which is in error by 2 in the 8 significant digit. 



Note that the ratio for the scaling limit of SAWs in a 10 x 1 rectangle is some two orders of magnitude larger than 
the corresponding result for Brownian motion. This can be qualitatively understood as SAWs having a greater 
tendency to go in straight lines than do random walks. 



4.4. Pulling a polymer from a wall 



During the past decade, force has been used as a thermodynamic variable to understand molecular interactions 
and their role in the structure of bio-molecules [39, 24j|42]. By exerting a force in the picoNewton range, one 
aims to experimentally study and characterize the elastic, mechanical, structural and functional properties of 
bio-molecules [8 . 

In [57] SAWs were used to model the situation in which a polymer is attached to a surface and pulled from that 
surface by an applied force. The situation is shown in Figure [7j Interactions are introduced between neighbouring 
monomers on the lattice that are not adjacent along the chain. The pulling force is modeled by introducing an 
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energy proportional to the :r-componcnt of the end-to-end distance. One end of the polymer is attached to an 
impenetrable surface while the polymer is being pulled from the other end with a force acting along the x-axis. 
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FIGURE 7. An interacting self-avoiding walk on the square lattice with one end attached to a 
surface and subject to a pulling force at the other end. Each step of the walk connecting a pair 
of monomers is indicated by a thick solid line while interactions between non-bonded nearest 
neighbour monomers are indicated by jagged lines. 

Boltzmann weights uj = cxp(— e/k-^T) and u — exp(— F '/k^T) conjugate to the nearest neighbor interactions 
and force, respectively, were introduced, where e is the interaction energy, &b is Boltzmann 's constant, T the 
temperature and F the applied force. For simplicity, we set e = — 1 and fee = 1. The relevant finite- length 
partition functions are 

Z N {F,T) = YjJ m u x = ^jC{N,m,x)u m u x , ( 12 ) 

all walks m,x 

where C(N, m, x) is the number of interacting SAWs of length N having m nearest neighbor contacts and whose 
end-points are a distance x = xn — x apart. The partition functions of the constant force ensemble, Zjq{F,T), 
and constant distance ensemble, Zn{x, T) — J2 m C(N, m, x)uj m , are related by Zjy(F, T) = J2 x Zjy(x, T)u x . The 
free energies are evaluated from the partition functions 

G{x) = -TlogZ N {x,T) and G(F) = -Tlog Z N (F, T). (13) 

Here (x) — dG Qp^ and (F) — dG g ^ are the control parameters of the constant force and constant distance 
ensembles, respectively. 

All possible conformations of the SAW were enumerated. The challenge facing exact enumerations is to increase 
the chain length. Using direct counting algorithms the time required to enumerate all the configurations increases 
as fj, N , where /i is the connective constant of the lattice (/i ss 2.638 on the square lattice). So even with a rapid 
increase in computing power only a few more terms can be obtained each decade. In |27| the number of interacting 
SAWs was calculated using transfer matrix techniques [22 . Combined with parallel processing these algorithms 
allowed the enumerations to be extended to chain lengths up to 55 steps, roughly doubling the previously available 
enumerations. 

In Fig. |8^, we show the force-temperature phase diagram for flexible chains. At low temperature and force the 
polymer chain is in the collapsed state and as the temperature is increased (at fixed force) the polymer chain 
undergoes a phase transition to an extended state. The transition temperature as plotted in Fig. [8^1 was found (in 
the thermodynamic limit N — > oo) by studying the reduced free energy per monomer. The most notable feature 
of the phase-diagram is the re-entrant behaviour-that is to say, the initial increase of force with temperature, 
prior to a decrease. 

The positive slope dF c /dT at T — confirms the existence of re-entrance in the F — T phase-diagram. The 
authors pointed out that the value of the transition temperature obtained in the thermodynamic limit and the 
one obtained from the fluctuations in non-bonded nearest neighbors (which can be calculated exactly for finite 
N = 55) gives almost the same value (within error bars of ±0.01). The fluctuations are defined as \ — (m 2 )~(m) 2 , 
with the k'th moment given by 

(m fe ) = ^m k C{N,m,x)Lo m u x /Y^C{N,m,x)u m u x . 

m,x m,x 
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a 




T T 



FIGURE 8. a The force-temperature phase diagram for flexible chains, b The fluctuations \ m 
the number of contacts vs. temperature at fixed force F = 1.0 for various values of the chain 
length N. 

In the panels of Fig. [SJd the emergence of two peaks in the fluctuation curves with increasing N at fixed force 
F = 1.0 are shown. The twin-peaks reflect the fact that in the re-entrant region as one increases T (with F fixed) 
the polymer chain undergoes two phase transitions. The importance of powerful enumeration data is highlighted 
by the observation that the twin-peaks are not apparent for small values of N. Many more details and comparison 
with experiments are given in |27| - our purpose here is just to show the applicability of SAWs to this problem. 

5. The scaling limit and SLE 

5.1. The scaling limit 

One topic we have failed to adequately address is the question of the scaling limit of SAW. An intuitive grasp 
of this concept can perhaps be gained by looking at the first two figures in this article. In the first figure, the 
effect of the lattice is clear. In the second figure, there is no obvious lattice, and indeed no way to tell that this 
is not a continuous curve. We formalise this notion as follows: Consider a smooth (enough) closed domain f2, 
with an underlying square grid, with grid spacing 8 as shown in Figure |5.1| Denote by fls that portion of the 
grid contained in £1. Take two distinct points on the boundary of Q labelled a and b. Now take the nearest lattice 
vertex to a, and label it as, and similarly b$ is the label of the nearest lattice vertex to point b. Consider the 
set of SAWs cj(Qs) on the finite domain £1$ from as to bs- Recall that S > sets the scale of the grid. Now let 
|w| be the length of a walk us € oj(£Is), and weight the walk by x'"'. The reason for this is that the walks are 
of different lengths, making the uniform measure not particularly natural. (There is also a normalising factor, 
which for simplicity we ignore). 
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Figure 9. Discretisation of domain f2. 



As we let S — > we expect the behaviour of the walk to depend on the value of x. For x < x c it is possible to prove 
that cos goes to a straight line as S — > 0. (Strictly speaking it converges in distribution to a straight line, with 
fluctuations O(VS).) For x > x c it is expected that uj$ becomes (again, in distribution) space-filling as S — » 0. 
But at x — x c it is conjectured that uj$ becomes (in distribution) a random continuous curve, and is conformally 
invariant. This describes the scaling limit. If this conjecture is correct, a second, pivotal, conjecture by Lawler, 
Schramm and Werner |31j is that this random curve converges to SLE 8 / 3 from a to b in the domain Q. These two 
conjectures must be considered the principal open problems in the field. If they could be proved, the existence 
and value of the critical exponents, as predicted by conformal field theory for two-dimensional walks would be 
proved. 



5.2. Schramm Lowner Evolution 



For an approachable discussion of SLE K , see Chapter 15 of [16] . Here we give a very minimal outline. Let H 
denote the upper half-plane. Consider a path 7 starting at the boundary and finishing at an internal vertex. 
Then H\7 is the complement of this path, and is a slit upper half-plane. It follows from the Ricmann Mapping 
Theorem that it can be conformally mapped to the upper half-plane. Lowner [29J discovered that by specifying 
the map so that it approaches the identity at infinity, the conformal map so described (actually a family of maps, 
appropriate to each point on the curve) satisfies a simple differential equation, called the Lowner equation. The 
mapping can alternatively be defined by a real function. This observation led Schramm to apply the Lowner 
equation to a conformally invariant measure for planar curves. That is to say, the Lowner equation generates 
a set of conformal maps, driven by a continuous real-valued function. Scramm's profound insight was to use 
Brownian motion B t as the driving functiorj^J So let B t , t > be standard Brownian motion on K and let k be 
a real parameter. Then SLE K is the family of conformal maps g t : t > defined by the Lowner equation 

This is actually called chordal SLE K as it describes paths growing from the boundary and ending on the boundary. 
If K < 4 then the path is almost surely a simple curve, in the upper half plane. Larger values of k lead to more 
complicated behaviour. 



It is the only process compatible with both conformal invariance and the so-called domain Markov property. 
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6. Conclusion 



Hopefully this rather vague description will convey the flavour of this exciting and powerful development in 
studying not just two-dimensional SAWs, but a variety of other processes, such as percolation, the random 
cluster model, and the Ising model. We refer the reader to [I] for greater detail of both SLE and these appli- 
cations. Despite these remarkable advances, we still have no real idea how to obtain comparable results for the 
3-dimensional mode0 

In this article I have only scratched the surface of this topic. More details on the mathematical aspects can 
be found in |33| and the recent reviews [21|4j. More information on numerical aspects and some applications, 
particularly to the SAP subset can be found in the monograph [16 . Another approach to this problem that has 
not been discussed is to simplify the problem so that it can be solved (see [IB] Chap. 3). Unfortunately most 
such simplifications involve rendering the model Markovian, which removes a significant feature. 

As can be seen, there are many challenging open problems in the field that might capture the imagination of 
mathematical physicists, computer scientists and probabilists. Perhaps the major open problem is to prove that 
the scaling limit of SAWs is described by SLEg,/ 3 . Even proving the existence of critical exponents for two- or 
three-dimensional SAWs would be a significant advance. 
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